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We describe recent progress on generation of gauge configurations using the Highly Improved 
Staggered Quark (HISQ) action that was designed by the HPQCD/UKQCD collaboration. The 
HISQ action requires two levels of smearing with a reunitarization of the links before the second 
smearing. We describe how we deal with the occurrence of occasional large forces arising from 
the reunitarization step. The MILC collaboration is currently generating ensembles with approx- 
imate lattice spacings of 0.15, 0.12, 0.09, and 0.06 fm, with the strange and charm quark masses 
close to their physical values and the mass of the light quarks to/ set to 0.2m s . We present recent 
results for pion taste splittings, light hadron masses, the static potential, the r\ c dispersion relation 
and the topological susceptibility. 
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1. Introduction 

The MILC collaboration has been using the asqtad action ^ for more than a decade for its 
generation of 2+1 flavor dynamical quark configurations. Over forty ensembles covering six lattice 
spacings from a 0. 18 fm to 0.045 fm and a wide range of quark masses have enabled us to explore 
both the continuum and chiral limits 

Recently, the Highly Improved Staggered Quark (HISQ) ^ action was developed by the 
HPQCD/UKQCD collaboration. This action requires two levels of smearing with a reunitariza- 
tion of the links before the second smearing. Although it is more expensive to simulate than the 
asqtad action, its reduced taste symmetry breaking makes the HISQ action an attractive possibility 
for future calculations. It also allows us to include a dynamical charm quark with reasonable accu- 
racy at current lattice spacings. The goal of this initial work is to generate several ensembles with 
a fixed ratio of light to strange quark masses and to study the approach to the continuum limit. 

Large forces arise during molecular dynamics updating when the determinant of a link at the 
first level of smearing is close to vanishing [Qj. This problem is dealt with by a slight change in the 
guiding action, the effect of which can be eliminated via the accept/reject step at the end of each 
trajectory. 

The MILC collaboration is currently generating ensembles with lattice spacings of about 0.15, 
0. 12, 0.09, and 0.06 fm, with the strange and charm quark masses near their physical values and 
the mass of the light quarks mi set to 0.2m s . In the future, we plan to generate additional ensembles 
with lighter up and down quark masses. We present recent results on pion taste splittings, light 
hadron masses, the rj c dispersion relation and the topological susceptibility. 



2. Implementation of the HISQ action 

The HISQ action was developed to go beyond the asqtad action in its improvement of taste- 
symmetry. Links in the Dirac operator are replaced by 

U' = ^ 2 ^^iU (2.1) 

where smearing level 1 denoted &\ is Fat 7; reunitarization is denoted by ; and smearing level 2 
denoted J? 2 is asq smearing. 

To calculate the fermion force, we must take the derivative of the action with respect to the 
gauge link. With multilevel smearing, we must use the chain rule []|, ^]. For the fermion force, we 
have: 

dS f _ dS f dX dW dV 

~dU ~ ~dX dW ~dV dU ' 
where U are the fundamental gauge links, V represent fat links after level 1 smearing, W are the 
result of reunitarizing V, and X are the links after asq smearing. Conveniently, code for three of the 
derivatives is the same as for the asqtad action (although the arguments may be different): 

d li ?L n „ 

dx ' dW ' du ' K ' 

After the first level of smearing, the link V must be projected into U(3). This required some 
new code. We define: W = VQ~ 1 / 2 , with Q = V^V. The method to calculate the inverse square 
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Figure 1: Time history of the minimum value of det | V | (upper plot) and the maximum magnitude of the 
fermion force (lower plot) [Q]. Note how the large spikes in the fermion force are correlated with near 
vanishing of det |V|. The lattice spacing is ~ 0.13 fm. 

root, based on the Cayley-Hamilton theorem, had already been applied by Morningstar & Peardon 
[0], and Hasenfratz, Hoffman & Shaefer 

Q~ 1/2 =fo + fiQ + f2Q 2 . (2.4) 

We just need to find the eigenvalues of Q. 

A significant issue in implementing molecular dynamics based on this algorithm is that there 
can be large changes in the action resulting from large forces during the integration of the equations 
of motion [Q]. In Fig. [I], we see that there is a correlation between the large forces and low values 
of det | V | . In Fig. ||, we show the fermion force rescaled by I /At for runs with three different 
lattice spacings. (These two graphs are based on exploratory runs details of which can be found 
in Ref. [Q].) It appears that as we approach the continuum, there are fewer large forces. As the 
fundamental links U are varying more smoothly, we should not be surprised that V is less likely to 
have a small determinant. 

We implement an "eigenvalue filter" to tame the large forces. If the smallest eigenvalue of a 
smeared link V is smaller than a certain cutoff, we make the replacement <2 1//2 - > (Q + Sl)^ 1 ^ 2 , 
where 8 is typically set to 5 x 10~ 5 . Although this modifies the guiding Hamiltonian, the integra- 
tion algorithm is still reversible and area preserving. As long as we use the RHMC algorithm [g], 
the accept/reject step at the end of each MD trajectory ensures the correct equilibrium distribution. 

3. Testing the HISQ action 

For our initial test of the HISQ action, we are generating ensembles with mi = 0.2m s and 
lattice spacings varying from a w 0.15 to 0.06 fm. Table 1 shows the masses, lattice dimensions 
and acceptance rate on our tuning runs. In Table 2, we consider the production runs at the two 
intermediate lattice spacings. This range of lattice spacings is comparable to what we have with 
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Figure 2: Time history of maximum value of the fermion force rescaled by the inverse time step for ex- 
ploratory runs with three different lattice spacings. See Ref. Q for run parameters. 
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Table 1: Parameters for tuning runs: coupling, dynamical quark masses in lattice units, lattice dimensions, 
approximate lattice spacing, and acceptance rate. 
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size 
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6.00 
6.30 


0.0102 
0.0074 


0.0509 
0.0370 


0.635 
0.440 


24 3 x 64 
32 3 x 96 


0.12 
0.09 


0.74 
0.77 



Table 2: Same as for Table 1, except for production runs. 



the asqtad action [Q] and we compare results for a number of quantities. We find that the pion taste 
splittings are reduced by a factor of 2.5 to 3, and that results with the HISQ action are comparable 
to those using the aqstad action at a lattice spacing of about 2/3 that used for the HISQ action. 

4. Initial results 

We present some initial results for the taste splittings of the pion, the spectrum of particles 
with light quarks, the r\ c dispersion relation and the topological susceptibility. In most cases, the 
results are compared with prior work on asqtad configurations and it is found that the results with 
the HISQ action look closer to the continuum limit. (We do not have results for the r\ c dispersion 
relation with asqtad quarks.) 
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4.1 Pion taste splittings 

It is well known that for staggered quarks, taste symmetry breaking is an important lattice 
artifact. A prime motivation for HISQ quarks is to improve the taste symmetry. In Fig. ||, we 
compare the pion taste splittings for the aqstad [Q] and HISQ actions. We see that for a w 0.12 
fm, the taste splittings for the HISQ action are somewhat smaller than for asqtad with a w 0.09 fm. 
Similarly, results using the HISQ action with a rs 0.09 fm are similar to results using the asqtad 
action with a m 0.06 fm. The vertical line to the left of this semi-log plot marked x3 is a ruler 
showing a factor of 3. We see that using the HISQ action reduces the taste splitting by a factor of 
2.5 to 3 compared with asqtad. 
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Figure 3: Calculation of taste splitting using the asqtad and HISQ actions. We show {M\ — Mo)r\ vs. 
(Xy(a/ri) 2 , where Mq is the Goldstone taste state. The diagonal line on this log-log plot shows the slope of 
the expected dependence of the taste splitting. 

4.2 Spectrum of other light hadrons 

We also have some preliminary results for the masses of p, nucleon and Q. . In each case, 
we find that using the HISQ action results in masses with smaller lattice spacing dependence, and 
the results are comparable to asqtad results at smaller lattice spacings. (The poster presented at the 
conference did not show results for p due to lack of space.) In Figs. |]and ||, we show results for the 
p, nucleon and Q. . The improvement for the p seems quite substantial. For the nucleon, we find 
that the HISQ action result for a sa 0.15 fm is quite a bit below the asqtad result with a 0.12 fm. 
The HISQ result for a sa 0.09 fm appears comparable to that for the asqtad action for a k 0.06 fm. 
The magenta curve represents a chiral perturbation theory fit to the asqtad results. For the Q. we 
also see a substantial improvement of HISQ action results over asqtad action results with a w 0.12 
and 0.09 fm. 
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Figure 4: Rho mass in units of r\ vs. the pion mass squared. 
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Figure 5: Nucleon mass (left) and il mass (right) in units of r\ vs. the pion mass squared. 

4.3 HISQ action and the charm quark 

The HPQCD/UKQCD collaboration has developed a procedure for tuning the coefficient of 
the Naik term to reduce discretization errors up to &{(am c ) A ). At tree level, the correction to the 
Naik term expressed in terms of the bare charm quark mass is [Q] 

21, , 2 327 , n4 15607 . . 6 73697 . xS 

£ = (am c Y^ (am c y (am c )° (am c r . (4.1) 

40 v 1 1120 v 1 268800 v ' 3942400 v ; 

Note that this formula differs from Eq. (24) in Ref. [||] because it is expressed in terms of the bare 
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mass appearing in the action. In Fig. 6, we show the dispersion relation for the rj c meson. For the 
three lattice spacings displayed, —0.355 < £ < —0. 123. For a « 0.06 and mi = 0.2m s , £ « —0.063. 
We see that for a lattice spacing as coarse as 0.15 fm, the deviation from the continuum dispersion 
relation is as large as 7%; however, when the lattice spacing is reduced to 0.09 fm there is less than 
a 2% deviation from the continuum for the range of momentum shown here. 
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Figure 6: Speed of light from the rj c dispersion relation vs. lattice momentum squared. 



4.4 Topological Susceptibility 

In Fig. 7, we plot the topological susceptibility for both asqtad and HISQ action simulations. 
We find that the HISQ action point for a « 0. 12 fm falls close to the asqtad action curve for a « 0.09 
fm. Also note that the HISQ action point lies far to the left of the corresponding asqtad action 
a ?s 0.12 point (the third green square from the left) since the HISQ action taste symmetry is so 
much better than for the asqtad action, and the horizontal axis involves the taste singlet pion mass 
as shown by Billeter, DeTar and Osborn [10]. The black points and curve are from a continuum 
extrapolation of the fits to the results for the three asqtad lattice spacings. 



5. Outlook 

Based on the encouraging results of our preliminary scaling study with m/ = 0.2m s , the MILC 
collaboration is planning to generate additional ensembles with mi = 0.lm s and 0.04/?^ for lattice 
spacings between 0.06 and 0.15 fm. These ensembles will allow us to study both the chiral and 
continuum limits. We expect to be able to study a wide variety of physics with these ensembles, 
for example, light pseudoscalar masses and decay constants, light quark spectroscopy, topology, 
spectroscopy of particles with a heavy quark, heavy-light meson decay constants and semileptonic 
decays, and quark masses. 
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Figure 7: Topological susceptibility vs. squared mass of the taste singlet pion. The symbol size has been 
increased for the points coming from the asqtad action that correspond to m/ = Q.2m s so that they can more 
easily be compared with the HISQ action point. These points are either the second or third point from the 
left, depending on the lattice spacing. 
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